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5j^^ ^ in D'' — !■ K^K^Tv decays are reported. These parameters can be used to improve the 



*~j ' determination of the unitary triangle angle 7 in B~ -^ DK~ decays, where D is either a D'' or a 

■ - - _D0 meson decaying to the same final state, and also in studies of charm mixing. The measurements 

of the coherence factor and strong-phase difference are made using quantum-correlated, fully- 
reconstructed D"D'^ pairs produced in e^e~ collisions at the ^(3770) resonance. The measured 

values are RffO i^^ — 0.73 ± 0.08 and 5^^ — (8.3 ± 15.2)° for an unrestricted kinematic region 

and Rk*k = 1.00 ± 0.16 and 5^''^ = (26.5 ± 15.8)° for a region where the combined A'^tt* 
invariant mass is within 100 MeV/c^ of the K*{892) mass. These results indicate a significant 
level of coherence in the decay. In addition, isobar models are presented for the two decays, which 
show the dominance of the K*{892)^ resonance. The branching ratio B(D° -> KgK+Tv~)/B{D° -^ 
K%K~ -K^) is determined to be 0.592 ± 0.044 (stat.) ± 0.018 (syst.), which is more precise than 
previous measurements. 

PACS numbers: 13.25.Ft,12.15.Hh,14.40.Lb,11.80.Et 

I. INTRODUCTION There are two possible decays, 1?° -^ K'^K-Tr+ and 

D^ — >■ KgK^Tr~ (charge conjugation is implied through- 
Decays of the D° meson to the K'^K^tt^ final state out the paper unless otherwise specified), where the first 
occur via a singly-Cabibbo-suppressed (SCS) transition. 



is decay is favored with respect to the second. These twin 
decays are of interest for a variety of reasons. Firstly 
when produced through a B^ — > DK^ decay, where D 
is either a D^ or D'^, they can be exploited to measure 
the Cabibbo-Kobayashi-Maskawa (CKM) phase 7 [l|. 
Secondly a time-dependent study can be used to deter- 
mine charm mixing parameters and probe for charm CP 
violation 2\. Finally, examination of the resonant sub- 
structure of the decays allows for tests of SU(3) flavor 
symmetry in charm meson decays [3|. 

This paper describes several measurements performed 
with data collected by the CLEO experiment that are 
valuable for each of these purposes. These are a first de- 
termination of the coherence factor and average strong- 
phase difference for the two decays, a determination of 
the relative branching ratio of the decays that is more 
precise than previous measurements [J|, and the presen- 
tation of isobar models to describe each decay. The co- 
herence factor Rj^o j^^ and the average strong-phase dif- 
ference Sj-,^ ^ ^ for D^ -^ KgK^T:^ are measured us- 
ing quantum-correlated, fully-reconstructed D^D^ pairs 
produced in e+e" collisions at the -0(3770) resonance. 
Knowledge of these parameters improves the sensitivity 
of measurements of the unitary triangle angle 7 using B~ 
meson decays to DK^ where the D decays to KgK~TT~^ 



or KgK^TT~ . To probe for the effects of physics beyond 
the standard model, it will be necessary to over-constrain 
the CKM quark-mixing matrix. An accurate measure- 
ment of 7 in i? — > DK decays is of particular impor- 
tance as the extracted phase is robust against new physics 
contributions, thus providing a standard model reference 
against which other measurements can be compared. The 
sensitivity to 7 comes from the interference of & — >■ m and 
6 — > c quark transitions. The weak phase between these 
two transitions is —7. The amplitudes are related by 
AiB- -^ D°K-)/A{B- -^ D^K-) = rBe^^^''-''\ where 
rs has a value around 0.1 and is the absolute ampli- 
tude ratio, and 6b is the strong-phase difference. Sen- 
sitivity to 7 is obtained by comparing the four separate 
decay rates involving B~ and _B+, and the two possi- 
ble D final states. This method is directly analogous 
to the ADS method [0], originally proposed to exploit 
Cabibbo-favored (CF) and doubly-Cabibbo-suppressed 
(DCS) decays such as Kn, Kinnr, and Ktttt''^ rather 
than the SCS modes under consideration here. In the 
case of decays to multi-body flavor specific final states 
it is necessary to take into account that the amplitude 
ratio of D° -> KIK~tt+ and D° -> K^K+t:- and the 
strong-phase difference will vary over phase space. Rate 
equations for different final states are given by ^>1 
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Here, Akoktit^ ("^/^o ^^ , ^k-k) i^ ^^^^^ amplitude for D^ — > 
K^K^TT^ at a point in multi-body phase space de- 
scribed by the square of the invariant mass of the 
K^K and Kt: pairs, (jti^o;^, ^Ktt) respectively, and 



42 



^drrir 



.dmi 



j \-^K0K^7r±{lT^K'^K^"'-KTTJ\ "'"•'k%K"-"''K-k- 

Analogous equations can be written for the decays Z?" — >■ 
K^ Ti^ 'ii~ TT^ and D^ — >■ K~7t~^7t'^ and the coherence fac- 
tors and average strong-phase differences for these final 
states have already been measured [7] . One difference be- 
tween the K^ir^, 



the KcK^n^ mode is that the value r^"^ is approx- 



i^T^Tr^TT+Tr and K^tt^tt^ modes and 

imately ten times larger than the corresponding param- 
eters for the other decays. The larger value for r^^ 
occurs because both final states occur via SCS transi- 



I 

tions, while in the case of the K^n^, K^tt^tt^tt^ and 
A'^TT^TT^ decay modes one transition is CF and the other 
DCS. The larger value of rjy means that the per-event 
sensitivity to 7 from the D" — >■ KgK~TT'^ decays [Eq. ([1])] 
is increased in comparison to D'^ — > K~it'^ decays, while 
the per-event sensitivity to 7 from the D^ — > K^K^ti~ 
decays [Eq. ([2])] is reduced in comparison to D" — >■ K^'k~ 
decays. The overall sensitivity to 7 remains significant 
if the coherence is large. Therefore knowledge of the co- 
herence factor in D^ — >■ K^K^tt^ decays is valuable, 
as it enables these channels to be added to the set of 
modes with which 7 can be measured using _B~ — > DK~ 
decays. The coherence factor and average strong-phase 
difference are also valuable for time-dependent charm- 
mixing measurements as described in [^. 



In a region around a resonance that is prominent in 
both decays, the coherence is expected to be close to one. 
A measurement in a bin of the Dahtz plot around such a 
resonance is therefore also useful as the higher coherence 
improves the per-event sensitivity in measurements to 
determine 7. In the case of D^ — > KgK^n^ decays 
a prominent intermediate resonance is the K*{892)^. 
Therefore, in this paper the measurement of the coher- 
ence factor and the average strong-phase difference is 
performed for both an unrestricted kinematic region and 
also for events where the K^tt invariant mass is within 
100 MeV/c^ of the nominal K*{892)^ mass [4]. The 
value of 100 MeV/c^ corresponds approximately to twice 
the natural width of the i^*(892)='= meson. 

In addition to the measurement of the coherence factor 
and average strong-phase difference, two further stud- 
ies are also presented. First, there is the determination 
of the isobar models for the two decays. These can be 
used to test the predictions of SU(3)-flavor symmetry 
in charm meson decays, in particular the relative am- 
plitud es an d phases of the modes Z?" — >■ K*^K^ and 
D" — >■ K*'^K'^ ji5(]. The isobar models are used to estimate 
systematic uncertainties for the other measurements pre- 
sented here and will also be of use as an input to sim- 
ulation when performing future measurements that ex- 
ploit Z?" — > KgK^n^ decays. Secondly, the ratio of the 
branching fractions for the suppressed with respect to the 
favored decay, Bko^t.^ B{D° -^ K%K+it-)/B{D^ -^ 
K^K~t:~^), is measured. The current world average is 
0.79 ± 0.19 [^], and a more accurate determination is 
possible using CLEO data. This is valuable as a better 
measurement of B^" k-w allows for a corresponding im- 
provement in the knowledge of r^^ , which is an im- 
portant input to the coherence factor analysis, and also 
to measurements of 7 and studies of charm mixing. The 
isobar model studies are based on data collected both 
at ^/s — 3770 MeV and from charm mesons produced 
in the continuum at higher energies; the branching ratio 
measurement exploits only the latter sample. 

This paper is organized as follows. Section |ll] gives 
an overview of the detector configurations used and the 
data samples that are used for all the studies presented in 
this paper. Section Hill describes the isobar models and 
Sec. IIVI details the measurement of the ratio of branching 
fractions. Section |V] describes the measurement of the 
coherence factor and the strong-phase difference; and in 
addition the sensitivity to 7 using these results is dis- 
cussed. A summary is given in Sec. IVIl 



II. DATA SAMPLES 

Two different data samples collected by the CLEO III 
and CLEO-c detectors are used for the studies in this 
paper. A total of 15.3 fb~^ of data was collected by the 
CLEO III detector. For these data the Cornell electron 
storage ring (CESR) accelerator operated at a range of 



energies from -y/i = 7.0 GeV to y/s = 11.2 GeV. For the 
data collected by the CLEO-c detector analyzed in this 
paper, the CESR energy was y/s — 3770 MeV and an 
integrated luminosity of 0.818 fb^^ of data was accumu- 
lated. The CLEO III/CLEO-c detector was a solenoidal 
detector which included a gaseous tracking system for 
the measurement of charged particle momenta and ion- 
ization energy loss, a ring-imaging Cherenkov (RICH) 
detector to aid in particle identification, and a Csl crys- 
tal calorimeter to measure the energy of electromagnetic 
showers. The tracking system was particularly important 
in defining the candidate signal sample for the measure- 
ments in this paper. The CLEO III tracking system Q 
consisted of a silicon strip vertex detector and a large 
drift chamber. The silicon was replaced by a six-layer 
wire vertex detector in the CLEO-c configuration [9|. 
The trackers achieved charged particle momentum reso- 
lution of 0.35% (0.6%) at 1 GeV/c in the 1.5 T (1.0 T) 
axial magnetic field used in CLEO III (CLEO-c) detec- 
tor, respectively. The detectors are described in detail 
elsewhere jlO| . 

The intermediate resonances contributing to the D" — > 
KgK^TT^ and D'^ — > K^K^tt^ decays are not identical, 
so in constructing isobar models it is necessary to know 
the flavor of the decaying I?°. The isobar models are de- 
termined from flavor-tagged data from the CLEO III and 
CLEO-c detectors. The data collected by CLEO III are 
flavor-tagged using the decay £)*+ — > £'°7r+, in which 
the flavor of the decaying neutral D meson is inferred 
from the charge of the pion. In the case of the data 
collected by the CLEO-c detector the interaction of in- 
terest is e+e~ — >■ V(3770) — >■ D°D° where one D me- 
son, the signal side, decays to K^K^tt^ and the other, 
the tag side, decays to either K^n^, K^n^n'^n^ , or 
K^TT^TT^. Reconstruction of both D mesons is referred 
to as a double tag. The flavor of the signal side D meson 
is inferred from the charge of the kaon on the tag side, 
taking advantage of the fact that the CF transition, for 
example D" — >■ K~Tr~^, is dominant. Corrections to this 
assumption from DCS decays are considered in the as- 
signment of systematic uncertainties. The measurement 
of Bj^o j^^ is performed using only the flavor-tagged data 
from 6lE0 III. 

The coherence factor measurement uses data only from 
CLEO-c. To maximize the sensitivity to the coherence 
factor a range of double tags are used. These include 
the flavor tags used for the isobar model and also dou- 
ble tags where the tag side D meson decays to a CP 
eigenstate or Kg ^tt+tt", which are states of mixed CP. 
Table [l] lists all the reconstructed double tags used for 
the coherence factor measurement. The neutral mesons 
in the final states are reconstructed via the following de- 



77, Kg — >■ 7r"^7r , u! — >■ 7r"^7r tt", 77 



77, 



The Kl meson 



cays: tt" ^ , ,, ^^g 

rj — > TT+TT^TT^, and 77' -^ ri{'yj)Tr^TT 

is not reconstructed directly. Its presence is inferred via 

a missing-mass technique as described in Sec. IIIBI The 

rest of this section describes the selection requirements 

placed on all the different data types used to perform the 



measurements presented in this paper. 



TABLE I. A list of the CLEO-c double tags that are selected, 
split by the category of the tag side D meson. All the listed 
double tags are used in the coherence factor measurement, 
but only the flavor double tags are used in the isobar model 
analysis. 

Tag group Tag side decays 



Flavor 
CP even 
CP odd 



K+K-,^+n-,Kl^\Klr^, 

Kluj, Kir,', K%^\° 

Kln°, Kin, Kico, KW, KW'n" 



Mixed CP K^TT+n-, KItt+tv- 



A. Data from CLEO III 

For the flavor-tagged D* data from CLEO III, various 
requirements are made on properties of charged tracks, 
composite particles and mass or vertex fits in order to 
achieve samples of D" — > K^K^tt^ virith high purity. 
The kaon and pion from the D meson and the tagging 
pion from the D* are required to pass close to the in- 
teraction point: within 5 cm along the beam axis, and 5 
mm perpendicular to it. The cosine of the angle between 
the track and the beam pipe is required to lie between 
—0.9 and 0.9, to ensure that the track lies within the 
detector acceptance. The momentum of the slow pion 
track is required to be between 0.15 and 0.50 GeV/c. All 
other tracks have a minimum momentum requirement 
of 0.2 GeV/c and a maximum momentum of 5.0 GeV/c. 
Charged pions are used to reconstruct the Kg candidates 
and the flight significance of the Kg candidate is required 
to be greater than 8, where the flight significance is given 
by the measured flight distance divided by the measured 
flight distance uncertainty. Kaons and pions are distin- 
guished by using information from the RICH detector if 
the momentum of the kaon track is above 500 MeV/c, 
otherwise ionization energy loss from the drift chamber 
is used. Kaon and pion candidates are combined with a 
Kg candidate to form a 1?° candidate. The slow pion 
candidate is combined with the D'^ candidate to form a 
D* candidate. In order to suppress combinatoric back- 
grounds and backgrounds from B meson decays, the mo- 
mentum of the D* is required to be greater than half 
its maximum possible value (for a given beam energy). 
Loose selection requirements are placed on the D'^ and 
D* fitted vertex x^ to remove combinatoric candidates. 

Signal candidates are chosen using the reconstructed 
D'^ meson mass rrij^o and Aitid — mu* ~ mjjo , where 
mo' is the reconstructed D* meson mass. The sig- 
nal region encloses a region 15 MeV/c^ either side of 
the nominal D^ meson mass, and 1 MeV/c^ either side 
of the nominal mass difference, where the nominal D'^ 
and D* masses are taken from [4]. These ranges are 



picked as they correspond approximately to a three stan- 
dard deviation window in either the Z?° mass or Amu 
projection, where the width of the signal is determined 
by a fit to each distribution. In each case, when the 
fit is performed on one variable the signal box cut has 
been applied to the other variable. The signal m^jo 
distribution is described by a bifurcated Gaussian and 
the background is described as a first-order polynomial. 
In AmD the signal is also described by a bifurcated 
Gaussian. The background shape is parametrized by 



(Amz,/a)2[l- 



3 — (Am_D— d)/c 



+b{AmD/d—l)] where a, b, c. 



and d are parameters determined by the fit. A correction 
is required for D^ — >■ K^K'tt'^ decays as Monte Carlo 
simulation studies indicate the presence of a peaking 
background. This background is the D^ -^ A'^tt+tt+tt^ 
decay misidentified as K'gK~T:^ and contributes approx- 
imately 2% to the peak. Monte Carlo simulation studies 
indicate no other significant peaking backgrounds for ei- 
ther decay. The signal yields from the fits are corrected 
to remove the D° — > X^tt+tt+tt^ contribution. The pu- 
rities of the samples in the signal region are evaluated by 
averaging the purity determined from the fits in both dis- 
tributions, which are consistent with each other. The pu- 
rity of the D" -^ KIK-t:+ sample is (92.4 ± 0.5)% with 
a signal yield of 640 ± 27 in the signal region after back- 
ground subtraction. The purity of the Z3" — > K^K^tt' 
sample is (90.0 ± 0.7)% with a signal yield of 406 ± 24 in 
the signal region. The given uncertainties on the purities 
and signal yields are statistical only. Figure [T] shows the 
D^ mass distribution and fit projections for favored and 
supressed decays after the signal window in Amu has 
been applied. The signal yield determination is also per- 
formed in the restricted kinematic region where the signal 
K^ and tt^ have an invariant mass within 100 MeV/c^ 
of the nominal iiT* (892)^ resonance. In this restricted 
region the signal yield is 445 ± 22 for D° -^ K^K-n+ 
events and 166 ± 14 for 13° — > KgK'^Tr~ events, and the 
purities are (93.8 ± 0.9)% and (90.9 ± 1.2)% respectively. 



B. Data from CLEO-c 

A variety of double tags, as listed in Table U are re- 
constructed in the CLEO-c data. The standard CLEO-c 
selection criteria for 7r+ 7r° and A'+ mesons are adopted 
and are described in |ll| . In addition to these criteria, 
the Kg vertex is required to be separated from the in- 
teraction point by at least two standard deviations and 
candidate Kg — >■ 7r"'"7r~ decays are required to have a 
mass within 7.5 MeV/c^ of the nominal mass. This cri- 
terion is reduced to within 5.0 MeV/c^ of the nominal 
mass when the Kg forms part of the KgK^n^ signal 
tag. This tighter cut is to reduce background from the 
£)° — > A'^TT+TT+TT^ decay faking the D° -> K'^K^tt^ 
signal. The 77 -^ 77 candidates are reconstructed in a 
similar fashion to tt*' candidates with the requirement 
that the invariant mass is within 42 MeV/c^ of the nom- 
inal mass; the same requirement is applied to ry — >■ tttttt^ 
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FIG. 1. The reconstructed mass distributions (points) of (a) D° -^ K%K-ti+ and (b) D° -> K%K+-k- in CLEO III data. The 
signal window selection in Amo has been applied. The total fit projection (dark line) and background (light line) component 
are also shown. The signal window in the reconstructed mass variable is indicated by the vertical lines. 



candidates. Candidates for uj — > tttttt^ decays are re- 
quired to be within 20 MeV/c^ of the nominal w mass. 
The 77' — >■ rjTTTT candidates are required to have an invari- 
ant mass in the range 950 to 964 MeV/c^. All nominal 
masses are taken from [^. 

For final states that do not contain a K^, two kine- 



matic variables are considered: 
candidate mass nibc = \/s/{4:c'^) — p^/ 



the beam-constrained 
where pd is 



the D candidate momentum, and AE — Ed — Vs/2, 
where Ed is the sum of the energies of the daughters that 
comprise the D meson candidate. Decays with correctly 
identified signal and tag decays peak at the nominal D^ 
mass in rjibc and zero in AE. Mode-dependent require- 
ments are placed on the signal and tag candidates such 
that AE is less than three standard deviations from zero. 
The double-tag yield is determined by counting events in 
signal and sideband regions of the [mbc{Dsig), mbc{Dtag)] 
plane. Figure [2] shows all the fully-reconstructed double 
tags on the [r7ibc(-Dsig),m.6c(-Dtag)] plane with the signal 
region marked S, and the sidebands illustrated. The rect- 
angular sidebands are populated by decays where only 
one of the D mesons has been correctly identified. The 
triangular sidebands are sparsely populated by combi- 
natoric background as neither D meson candidate has a 
mass consistent with a D meson. The remaining narrow 
diagonal sideband is populated with two types of events. 
One of these is continuum background where the orig- 
inal interaction was e^e~ — >■ qq, where q =u, d, or s. 
The other type of background present in this sideband is 
where tracks have been assigned to the wrong D meson. 
Also shown in Fig. [2] are the projections of the signal D 
meson mass in different tag groups. 

The different sidebands contain contributions from dis- 
tinct types of combinatorial background. The yields 



in these sidebands are scaled and subtracted from the 
yield in the signal region. There are some backgrounds 
that peak in the signal region and hence cannot be 
determined from the sidebands. These peaking back- 
ground contributions are estimated from Monte Carlo 
samples. EVTGEN [l2] is used to generate the de- 
cays and GEANT [3 is used to simulate the CLEO-c 
detector response. A generic sample of Monte-Carlo- 
simulated data of -0(3770) — >■ DD decays, with a size 
corresponding to an integrated luminosity twenty times 
that of data is used. The subsequent decay of the D" 
or D^ meson is to all known decay channels. A further 
sample of simulated continuum data of e+e" — )■ qq where 
q is either the u, d, or s quark, with a size corresponding 
to an integrated luminosity five times that of the data is 
also analyzed. All signal events, determined from truth 
information, are removed, and all background events that 
pass the selection criteria for any double tag are used to 
determine the peaking background for that double tag. 
The background contribution from the continuum is de- 
termined to be small in comparison to that from generic 
DD decays. The samples are combined by scaling the 
small continuum contribution to the integrated luminos- 
ity size of the generic sample. The statistical uncertainty 
due to the Monte Carlo background yields is very small 
in comparison to the uncertainty on the data yields. The 
uncertainty on the scale factor used to determine the 
peaking background yields from Monte Carlo simulation 
is 10%. This is determined by comparing all the double 
tag sideband yields in data to those in Monte Carlo to 
account for differences between Monte Carlo simulation 
and data. 

The largest peaking background to D'^ —> K^K^tt^ 
decays that remains after selection is from D'^ — >■ 
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FIG. 2. CLEO-c data distribution of fully-reconstructed double tags on the [mbc{Dsig),mbc{Dt!!,g)] plane (a) where the signal 
region is labeled S and the other indicated boxes are the sideband regions. The projections of the mtc{Dsig) are shown separately 
for candidates that are tagged by (b) CP eigenstates, (c) the flavor states, and (d) Kg-K'^Tv~ . 



K~TT^Ti^Ti^ decays where a tt+tt" pair forms a K'g can- 
didate. The background-subtracted signal region yields 
and purities from flavor tags are given in Table |ll] where 
the errors include the statistical uncertainties from the 
signal, sideband and Monte Carlo yields. Yields of other 
fully-reconstructed CP double tags and the K^rni tag 
are given in Table Hill For the X^tt+tt^ tag, the yields 
are required separately for 16 specific regions on the 
i^T^vr+TT" Dalitz plot. These regions are the equal phase- 
difference /S.5d bins reported in |14|. which are defined 
using the isobar model reported in [15| . Signal yields in 
these bins are determined and range between 0.7 ± 1.0 
and 7.8± 2.8. These uncertainties include the statisti- 
cal uncertainties on the signal region yield and the sub- 
tracted background. 

For final states that include a K^ meson, the missing 
mass squared recoiling against the fully-reconstructed D 
candidates and the particles from the other D decay con- 
taining the K^ meson is computed. Events consistent 
with the mass of the K^ meson squared are selected. 



TABLE II. Background-subtracted yields and purities of 
D" -^ KlK'^-K^ for each of the flavor tags in the CLEO- 
c data. The uncertainties include the statistical uncertainties 
of the yields in the signal region, sideband regions and Monte 
Carlo yields that determine the peaking backgrounds. 



Tag 



Yield 



VgK TV ' 

Purity % 



Yield Purity% 



K-K 122.0±11.5 93.2±0.5 
K-K-K-K 161.8±13.4 89.8±0.5 
KTi-n° 221.4±15.6 91.0±0.4 



80.6±9.0 99.5±0.1 
101.5±10.2 97.6±0.3 
120.8±11.1 97.4±0.2 



This technique was introduced in [16|. Further selec- 
tion requirements are as those listed in |14| . The com- 
binatoric background yield in the signal region is esti- 
mated from the population in the lower (L) and upper 
(H) missing mass squared sidebands. Information from 
the generic simulation is used to determine the relative 
composition of the sidebands and the signal region to es- 



TABLE III. Background-subtracted yields for fully-reconstructed CP tags and KgTVTT tag over all Dalitz plot bins. The 
uncertainties include the statistical uncertainties of the yields in the signal region, sideband regions and Monte Carlo yields 
that determine the peaking backgrounds. 



KK 



K^sVill) K'^rijnnn'') KW K^s^^ 



Yield 1.7±2.5 3.2±3.1 10.6±3.5 22.8±4.9 12.3±3.8 6.9±2.( 



0.4±1.0 0.8±1.0 60.7±8.1 



timate better the combinatorial background. There are 
also peaking backgrounds that only appear in the sig- 
nal region of the missing mass squared distribution. The 
largest contribution of peaking background arises from 
Z?" — >■ KgX decays, where the Kg is not reconstructed 
and the decay is mistaken for D^ -^ K^X. As an ex- 
ample. Fig. [3] shows the missing mass squared spectrum 
for the K^TT^TT~ tag in data and the background Monte 
Carlo distribution scaled to the integrated luminosity of 
the data. In the signal region the peaking contribution 
of the i^r^TT+TT" decays can be seen in the background 
Monte Carlo. The lower and upper sidebands are chosen 
to avoid any peaking structures that do not contribute to 
the background under the peak such as the contribution 
of D^ -^ K^ir^TT^ that is observed to the left of the lower 
mass sideband in Fig. |31 



tion and reconstruction efficiency for the K'^rj^TT^TT^Tr'^) 
tag, the Klri{TT^TT~TT'^) tag is not worth pursuing. Al- 
though several of the yields are close to zero the corre- 
sponding tag categories are retained in the analysis as 
the sensitivity to the coherence factor comes from the 
combination of many tags. The selection is repeated for 
all tags where the invariant mass of the Kg and tt from 
the signal decay is required to lie within 100 MeV/c^ of 
the K*{892)^ nominal mass. The yields in this restricted 
region typically lie between 40% and 80% of the reported 
unrestricted yields, depending on tag type. The selected 
samples in the restricted region have similar purities to 
those in the unrestricted region. The efficiency-corrected 
yields of these samples are given in Sec. IV Dl 



III. ISOBAR MODELS 
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The intermediate resonance structure in the decays 



D° -^ K^^K-jT+and Z?" 



KgK^TT can be visualized 



through the distribution of the data over the Dalitz plot. 
The Dalitz-plot distributions for each flavor-tagged data 
set are shown in Fig. U where all three flavor tags used 
in CLEO-c are plotted together. The horizontal band in 
each plot is a clear indication of the presence of an in- 
termediate K*{892)^ resonance. The remainder of this 
section describes the construction of isobar models for 
the favored and suppressed modes that parametrize the 
relative amplitudes and phases of all intermediate reso- 
nances that contribute to the decay. 



FIG. 3. The distribution for the missing mass squared for the 
KItv'^-k" tagged candidates. The points show the data and 
the solid line is the Monte Carlo background prediction, scaled 
to the data luminosity. The region between the two straight 
solid lines is the signal region, marked S, and the regions 
between the two dotted lines are the sidebands, marked L 
and H, used to determine the combinatorial background. 



The background-subtracted signal yields of the CP 
tags with a K^ meson in the tag side and the K^tttt tag 
are given in Table ITVl For the K^tttt tag, the candidate 
events are separated into regions on the Dalitz plot and 
the background-subtracted yields in each region range 
from 0.8 ± 1.4 to 24.8 ± 5.7. The purity of X^tt" tagged 
candidates is 67%. Given the very low yield for the can- 
didates tagged by ^^^'7(77) ^^'^ ^^^ lower branching frac- 



A. Signal model for the amplitude fit 

The signal isobar model considers the decay D^ — > 
CR, R — > AB, where A, B, C are pseudoscalar final 
state particles, and R is an intermediate resonance with 
spin 0, 1, or 2. The corresponding four-momenta are po 
for the initial Z?, pn for the resonance, and pa,Pb,Pc for 
the final state particles. For the invariant mass of the res- 
onance (which, due to its width, will vary event by event), 
the symbol itiab is used, with m^g = Pr = (PA + Pb)^- 
For the "nominal" mass of the resonance the symbol nifj 
is used. The matrix element for a D" — >■ CR, R -> AB 
decay chain is of the general form for resonance Ri 



M,=F, 



n, 



D- 



rrir 



m 



AB 



imjiT 



-Fn 



AB 



(5) 



TABLE IV. Background-subtracted signal yields for the CP tags with a K1 and the K^tvtt tag over all Dalitz plot bins. The 
uncertainties include the statistical uncertainties of the signal region and subtracted background estimates. 



Klu Ktrii-f-y) 



"kW 



Kt TV TT 






Yield 22.0±5.8 2.9±2.7 -0.1±1.8 0.3±1.0 2.8±2.0 150.9±14.2 
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FIG. 4. Dalitz plot distributions for (a)[(c)] D" -^ K%K~-r+ and (b)[(d)] _D" ^ _ft:gii'+7r- for flavor-tagged CLEO III [CLEO- 
data. 



where Fj:, and F^ are orbital angular momentum barrier 
penetration factors, l/{mJi—m\g — imiiTAB) is the line- 
shape of the resonance (in this example Breit-Wigner; 
other choices are possible), and J7j is the spin factor for 
a resonance of spin J. The running width. Tab is given 
by 



Tai 



f PAB 
\ Pr 



2,7+1 



rriAB 



Fl 



(6) 



where V^ is the natural width, pab is the momentum 
of one of the daughters in the AB rest frame, and Pr is 
the momentum of one of the daughters in the resonance 
rest frame. The spin factors are calculated as described 
in IITI, and are 



"2 - {p'd+P^'Wd +Pc)T^.ucp{v%-p'l){p% 



-p^)/(GeV)^ 



(7) 
(8) 
(9) 



r 



TABLE V. Particle ordering in spin 1 resonances. 



Decay 


A 


B 


C 


D" -^ K*+K-, K*+ -^ Kl;TT+ 
D° ^ p-TT+, p- ^ K-K°s 


K- 


K- 


K- 

7r+ 



TABLE VL The Blatt-Weisskopf penetration factors for the 
mesons of different spin. The momentum of either daugh- 
ter in the meson rest frame, Pr and the momentum of either 
daughter in the candidate rest frame pab, differ due to the 
use of the mass of the meson or the invariant mass from the 
two-daughter combination. 



Spin 


Fr 





1 


1 


1 i+tc^pP 


2 \ 


1 9+3TC2pp+-R4p4 


/ 9+3-R^pj^B+TZ'ip\g 



where T^^ = -g^" + Pt^Pa/m% and T^i.q/3 = 
(1/2) {T^c.T,0 + T^pT,^) - (l/3)r^,r„^. The sign of the 
spin-1 factor depends on the ordering of the resonance's 
decay products. The ordering choice for the different 
spin factors is given in Table |Vl Care is taken to en- 
sure that CP-conjugate amplitudes such as 13 — >■ K*K^ 
and D — >■ K*K^ are assigned the same spin factor (with 
the same sign). The Blatt-Weisskopf penetration fac 



tors, Ff> and Fu are described in [18| and are listed in 
Table IVII They rely on one free parameter 7?,, the me- 
son radius and on the momentum of the decay particles 
in the rest frame of the parent. For consistency with 
other experiments (la| TZ is fixed to be 5 (GeV/c)""'^, and 
1.5 (GeV/c)^^ for the D^ meson and intermediate reso- 
nances, respectively. 

The signal model is made from a sum of intermediate 
resonances so that the total matrix element is given by 



7Wtot = ^%e''^^A^j, 



(10) 



where Mj is as defined in Eq. ([5]), and aj and <t>i are 



the relative amplitudes and phases between different res- 
onances. There are a number of intermediate resonances 
which could contribute to the decays D" — >■ K^K^ir^and 
D° -^ K%K+-K-. The presence of the K*{SQ2)^ reso- 
nance is clearly visible in the Dalitz plot distributions 
for both the favored and suppressed decays. For this 
resonance a,- is set to one and dj to zero, as all other 



amplitudes and phases can be measured relative to this 
component. The Breit-Wigner lineshape is used for this 
resonance. The presence of the K^{\A?>Q)'^ is also re- 
quired. It is described with the LASS lineshape [2^ with 
the modifications described in [15|. This parametriza- 
tion is used as it includes an S-wave contribution and 
can account for the broad shape of any non-resonant 
component. The finite sample size limits the number of 
contributions that can be fitted to the data, without a 
large number of components having fit fractions consis- 
tent with zero, where the fit fraction Fi for resonance i 
is defined as 



F^ 



I '^^K"K'^^Kj^^imKOK'^K^)\^' 



I dmj^o^dmj^^ \ Y2=i ^k{ 



2 2 



(11) 



KttI 



The number of intermediate resonances fit to the favored 
decay channel is six as this usually produces a fit where 
no more than two resonances have fit fractions consistent 
with zero. The number of intermediate resonances fit to 
the suppressed decay channel is five, due to the lower sig- 
nal yield of this decay. Hence in addition to the if* (892)^ 
and iiTp (1430)* a further four (three) resonances are 
chosen to make the favored (suppressed) model. These 
resonances are taken from all known possibilities that 
decay into either if^TT, if^Tr^, or K^K^. These are 
the ao(980,1450)+, a2(1320)+, K*(892, 1410, 1680)0/±, 
is:^(1430)0/±, iiTs* (1430)0/±, and p(1450, 1700)+, where 
for each resonance the mass and width are taken from 
the PDG Q- The Breit-Wigner lineshape is used for all 
these resonances except for the iirQ(1430)+ which uses 
the lineshape in [l^, and ao(980) where the Flatte [21| 
lineshape is used. 



B. Detector efHciency and background modeling 

The detector efficiency, which varies over the Dalitz 
plot, and the distributions of the background events need 
to be taken into account when fitting the data to ex- 
tract the amplitudes and phases of each intermediate 
resonance. The models for the detector efficiency are 
determined from Monte Carlo simulation. The efficiency 
shapes for CLEG 111 and CLEO-c data differ due to the 
different D'^ momentum distributions and detector con- 
figurations. The simulated events are generated non- 
resonantly and hence are distributed evenly over the 
Dalitz plot. The shape of the relative efficiency is largely 
fiat although a reduction in efficiency is seen in regions 
where one of the particles has low momentum. The sim- 
ulation samples for CLEO-c are large. The probability 
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bands are not used due to insufficient yields to describe 
the shape on the Dahtz plot. 



C. Dalitz plot probability density function 

The data are fit using a similar method to previous 
CLEO analyses [12]. In order to describe the event den- 
sity distribution on the Dalitz plot a probability density 
function is used. For the isobar model fit the PDF is 
given by 



(1 - fsis,j)Bj{mj^oj^,mJi^), 



(12) 



density function (PDF) can be integrated using the sim- 
ulated events, and as they are distributed according to 
the efficiency shape there is no need to parametrize the 
CLEO-c efficiency. The Monte Carlo statistics available 
for CLEO III are not as large and in this case the effi- 
ciency shape is parametrized by a second-order polyno- 
mial in ("1^0^, m^^). 

The flavor-tagged data samples have a purity of 90% 
or greater as described in Sec. IIIAI and Sec. IIIBI For 
CLEO III data the background shape is determined from 
three sidebands in data. Two of the sidebands require 
Amu to be less than 1 MeV/c^ from its nominal value 
and either 1.7 < 171^0 < 1.8 GeV/c^ or 1.9 < rrijyo < 1.93 
GeV/c^. The second sideband has a more restricted 
mass range in order to avoid contamination from D^ — >■ 
7r+7r~7r+7r~ events that peak at higher reconstructed D" ^ 2 \ / 2 2 \ 

mass. The third sideband requires m^o to be within 15 ^h*^^*^ ^i^K«K^ ^kJ and S(m^o;^, m^^) represent the 
MeV/c^ of its nominal value and 0.148 < Am^i < 0.155 
GeV/c^. The candidate events from all three sidebands 
are combined as there are insufficient data for separate 
fits. The parametrization of sideband data over the 
Dalitz plot is determined empirically. It is modeled as 
a sum of four incoherent pseudo-resonances where the 
masses, widths and spins of the pseudo-resonances are 
adjusted until a good fit is achieved. Maximum likelihood 
fits are performed on background samples to determine 
the background shapes. The free parameters in these fits 
are the relative contributions of the pseudo-resonances. 
The background shapes are determined separately for the 
favored and suppressed decays as there are different con- 
tributions from misreconstructed D meson decays. For 
CLEO-c data the background shape is determined from 
the combined Monte Carlo samples described in Sec. Ill Bl 
in the same way as the CLEO III samples. The data side- 

I 



functions that describe the efficiency and background 
shapes, respectively, across the Dalitz plot, /sjg is the sig- 
nal fraction and j denotes either the CLEO III or CLEO- 
c data. The matrix element M is defined in Sec. IIII Al 
The PDF components e, |Mp, and B are all normalized 
individually over the Dalitz plot, which leads to V being 
normalized to unity as well. In the maximum likelihood 
fit to determine the isobar model the only free parame- 
ters are the real and imaginary part of each contributing 
intermediate resonance except the iiT* (892) where the 
amplitude and phase are fixed to one and zero, respec- 
tively. The background shapes and CLEO III efficiency 
shapes are fixed from fits to the individual relevant sam- 
ples described in Sec. IIIIBl and the CLEO-c efficiency is 
taken from the large Monte Carlo sample as described in 
Sec. IIIIBl The total likelihood to determine the isobar 
model is given by 



^- n 



n/sigjl-^("^KgK,»'"^l:7r,i)P£i(™XgJ<-,i'™|-7r,^) + (1 - /sigj)-Sj(m^0^,i,W^7r,*) 



j=CLEO-c, 
CLEO III 



(13) 



where i labels each event, in data sample j. 

Using the efficiency and background model PDFs, the 
likelihood given in Eq. (|13l) is constructed where the 
matrix element contains a combination of resonances. 
All possible isobar models, constructed from the full set 
of considered resonances, that have the allowed total 
number of components, and include the iiT* (892)^ and 
iirQ(1430)^ contributions, are tested when fitting to the 
data. 

The contributions of the resonances in the K^i:^ chan- 
nel are expected to be suppressed in comparison to reso- 
nances in the K^tt or KgK channel as this process pro- 
ceeds by a suppressed W boson exchange. Hence models 



I 

where there are large contributions in the K^tt channel 
are considered unphysical and are discarded. Also dis- 
carded are models where the total fit fraction is greater 
than 150% as these fits usually contain contributions that 
are highly correlated with each other and contribute very 
little independently to the fit. 

To estimate the quality of the fit a binned x^ quantity 
is computed. The Dalitz plot is split into square bins of 
size 0.18 GeV^/c"* vs. 0.18 GeVVc''. There are a total of 
67 bins which contain at least part of the allowed phase 
space. Bins that contain less than 10 events are merged 
with the neighboring bins until the content is greater 
than 10. The x^ is defined as (A'obs~-^cxp)^/-^cxp, where 
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A^obs and A^cxp are the number of observed and expected 
events, respectively. For a given fit, the ^ per degree 
of freedom (dof) is evaluated, and used to discriminate 
between different fit models. 

Models that contain the K*[%%2)^ resonance are inter- 
esting as they can be used to test SU(3)-flavor symmetry 
predictions [^ . As the best-fit models in both the favored 
and suppressed case do not contain this contribution, an 
additional model is constructed for each decay with this 
resonance required to be present. The models with the 
iiT* (892)'^ resonance are labeled as model 1 and the best- 
fit models are labeled as model 2. Tabic FVlIl summarizes 
the resonance contribution, the x^/dof, and the figures 
and tables that display the results for each model. 



D. Systematic uncertainties 

The systematic uncertainties for the measured ampli- 
tudes and phases are determined by performing fits with 
the same resonances as the chosen models in an alter- 
nate manner. The shift between the alternate and de- 
fault fit result is taken to be the systematic uncertainty 
on that parameter. The relative sizes of the systematic 
uncertainties are discussed in relation to the statistical 
uncertainty on the fitted amplitudes and phases and the 
fit fraction. 

The leading systematic uncertainty arises from the 
choice of values for the meson radii in the Blatt-Weiskopff 
penetration factors. In the alternate fit these are changed 
to 1.5 (GeV/c)-i and 0.5 (GeV/c)-i for the D° meson 
and intermediate resonances, respectively. It is observed 
that the changes in some of the fitted parameters in the 
favored models are up to three times the size of the statis- 
tical uncertainty, whereas in the suppressed models the 
changes are of the same order as the statistical uncer- 
tainty. This is understood as the favored models have a 
larger contribution from resonances where the pole mass 
lies outside the phase space. The change in the value of 
the penetration factor increases with increasing distance 
from the pole mass, and hence this effect is a larger source 
of systematic uncertainty for the favored model. 

The other leading systematic uncertainty is attributed 
to uncertainties in the measured masses and widths of the 
resonances, and the parameters of the iiTg (1430)^ line- 
shape. The masses and widths of the resonances, and the 
parameters of the iiTg (1430)^ lineshape are varied, each 
in turn by one standard deviation from their measured 
values. The shifts on the fitted amplitudes and phases 
are added in quadrature for each alteration. Overall the 
total systematic uncertainty arising from these changes 
is of order half the statistical uncertainty. 

The PDF used to arrive at the default models neglects 
the DCS amplitude that is present in the flavor tag de- 
cays. Although the DCS rate is much smaller than the 
CF rate, the interference between the two amplitudes is 
a significant effect in the quantum-correlated CLEO-c 
data. To determine the size of this effect the data are fit 



with a modified PDF for the CLEO-c data. In the case 
of the favored data, the alternate PDF uses a total am- 
plitude given by Aot = .4fav + 0.05 x 0.7 x Aupe*"^. The 
parameters of the suppressed model are fixed to their de- 
fault values, while the parameters of the favored model 
are free. The factor 0.05 arises from the approximate 
ratio of DCS to CF amplitudes in neutral D meson de- 
cays to K-K. There is, however, dilution as the Kttttt: 
and i^TTTr" tags are not fully coherent. The factor of 
0.7 accounts for the weighted average of the fully coher- 
ent Ktt decay and the coherence factors for the other 
two tags measured in [7[. Additionally there is an un- 
known absolute phase 6 between the favored and sup- 
pressed amplitude models. To account for this, the fit 
is repeated assuming a phase difference of zero, 7r/2, 
TT, and 37r/2. The largest shifts in the fitted param- 
eters are observed for a phase difference of zero and 
these are taken to be the systematic uncertainty. In 
general the systematic uncertainty due to this effect is 
about 25% of the statistical error. In the CLEO III 
data there is no quantum-correlation, however mistags 
occur when a track of the wrong charge is misidenti- 
fied as the pion from the D*. The mistag rate is de- 
termined by measuring the relative rate of, for exam- 
ple, Z?" — > K^TT^TT^TT— and D'^— > K^tt^tt^tt^ candi- 
dates associated with a D*'^ -^ D'^tt^ decay chain, and 
attributing to mistags the excess of doubly-Cabibbo- 
suppressed candidates that is observed above expecta- 
tion. The mistag rate is found to be around 0.6%. The 
systematic uncertainty due to this is determined to be 
negligible in comparison to other possible sources of bias. 

The systematic uncertainty due to the background 
fraction is evaluated by modifying the assigned value 
of the purity in the fit PDF by the uncertainty on 
its measured value. In addition, alternate background 
models are determined using a fit to the background 
samples with fewer pseudo-resonances. The alternate 
parametrizations are used in the fit to signal data, and 
the observed shifts are used as the systematic uncer- 
tainty. The uncertainty on the efficiency shape over the 
Dalitz plot is assessed using alternate parametrizations 
for the efficiency shape. In the case of CLEO III a third- 
order polynomial is used instead of the default second- 
order polynomial. In the case of CLEO-c the efficiency 
is parametrized by a fourth-order polynomial instead of 
using the Monte-Carlo-simulated sample directly. The 
uncertainties on efficiency, purity, and background shapes 
each contribute an uncertainty of approximately 10% of 
the statistical uncertainty. The default fit does not take 
into account the detector resolution on the Dalitz plot, 
and hence all data points are slightly shifted from their 
true positions. To account for the resulting uncertainty, 
the Kg, K, and tt momenta are smeared in data, accord- 
ing to parameters determined from Monte Carlo simu- 
lation, and the fit is performed on the smeared sample. 
This systematic uncertainty is negligible in comparison to 
other contributions. The validity of the fitter is assessed 
via fits to many samples of simulated data generated from 
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TABLE VII. A summary of the resonance contributions and x^/dof for each model. The figures showing the projection and 
the table listing the fit information are also given. 



Model 



Resonance 

A'*(892)+, A'*(892)", K'5{U3Q)+, K*{Umf, K*{1680)+, 



X /dof Projection Results 



Favored 1 

Favored 2 iC*(892)+, ^(1700)", is:o(1430)+ 

Suppressed 1 _K-*(892)-, _R'*(892)°, A'o*(1430)" 

Suppressed 2 j<"'(892r, A'2*(1430)", Kt^jlAJQ) 



ii'*(1410)°, A'*(1680)^ 

ao(1450)+, p(1700)+ 

, ao(1450)+, p(1700)^ 



ao(1450)- 


42.2/26 


Fig. [5] 


Table K] 


ao(1450)~ 


37.3/26 


Fig.E] 


Table IXII 




50.9/30 


Fig. [71 


Table Kill 




46.5/30 


Fig. El 


Table IXIIll 



the PDFs. A small uncertainty of order 10% of the statis- 
tical uncertainty is assigned due to small residual biases 
between the generated and fitted parameters. 

A summary of the systematic uncertainties for each ef- 
fect as a fraction of the statistical uncertainty on the am- 
plitudes, phases, and fit fractions for models 1 are given 
in Table Ivnll and HXl The relative sizes of the systematic 
uncertainties are similar for models 2. 



E. Results of the amplitude studies 



The results of the fit are listed in Tables IX HXilIl These 
tables show the amplitude, phase and fit fraction for each 
contributing resonance, and the statistical and system- 
atic uncertainty for each of these parameters. The fit 
projections for all these models are given in Figs. [5HS1 

The fit projections show good agreement with the data. 
The difference in fit projections between model 1 and and 
model 2 in either of the decay modes is very small. The 
iir*(892)^ component, clearly visible in the data, is seen 
to contribute significantly to all four models. This large, 
common resonant structure gives the expectation of sig- 
nificant coherence between the two decay modes, particu- 
larly in the region of the ii'*(892)^ resonance. The data 
for the suppressed decay show a broad extended shape 
either side of the 7^* (892)^ resonance in the ti^q^ pro- 
jection which is not present in the favored decay. This 
difference is apparent in the models as the favored models 
have a very small contribution from the i4r*(1430)° res- 
onance which is described by a broad S-wave lineshape, 
while in the suppressed models this is a significant con- 
tribution to the overall model. In the favored models the 
fit fraction of the K*{1A\0)^ resonance is large, which 
is unexpected as the K^tt^ channel proceeds by a sup- 
pressed W boson exchange. However the large fit fraction 
can be understood by considering the destructive inter- 
ference between K*{lAlQf and K*{IQSQ)^ which is also 
significant. This interference between two resonances can 
be quantified by calculating the interference fit fraction 
between any two resonances (j, /), and is given by 



l3,l = 



ReJdml.oj^dmJi^aje''''^Mja^e-^'^'Mi 



S dm\o^dm\^\YJl=iMk{ 



m' 



:.™k)l' 



(14) 



For these two resonances the interference fit fraction is 
-21% (-27%) for model 1 (2). The full tables of inter- 



ference fit fractions are provided in Appendix A. There is 
similarity between models 1 and 2 for each decay mode, 
as the contribution of the if* (892)° is small. From SU(3) 
flavor symmetry predictions 0] , the expected ratio of am- 
plitudes between D° -^ K*°K° and D° -^ K*+K~ is 
0.138±0.033. Using favored model 1 we observe a ratio 
of 0.2 ± 0.1. Although this is consistent with the predic- 
tion, it should be noted that the fit fraction of iir*(892)° 
is small and consistent with zero at the one standard de- 
viation level. An isobar model determined with a larger 
sample than analyzed here will be required to test the 
prediction with more meaningful precision. 



IV. MEASUREMENT OF THE RATIO OF 

BRANCHING FRACTIONS OF SUPPRESSED 

AND FAVORED DECAYS 

The measurement to determine the ratio of branching 
fractions uses only CLEO III data to avoid the complica- 
tions of accounting for the quantum correlations present 
in the CLEO-c data. The ratio ,B^o^^is given by 



B 



K°Ktt 



N{D° 



KlK^ 



e^a_ 



Oj^-, 



N{D° -> K'^K 



e_D"- 



(15) 



where N is the number of signal events for each decay, 
and e is the efficiency. 

The ratio of the efficiencies is not exactly unity as the 
distributions of the events over the Dalitz plot are dif- 
ferent in the favored and suppressed decays and the effi- 
ciency over the Dalitz plot is not uniform. The ratio of 
the efficiencies is found by taking the ratio of the inte- 
gral over all of phase space for the default isobar models 
multiplied by the efficiency model, for the favored and 
suppressed decay. The efficiency ratio is determined to 
be 0.940 using favored model 1 and suppressed model 
1. The uncertainty is assessed by altering the detector 
efficiency parametrization and using various alternative 
models arising from the fitting procedure that give a rea- 
sonable fit to data and leads to an uncertainty of ±0.029 
on the correction factor. Using the tagged CLEO III 
yields and the necessary corrections, Bj^oj^^ is measured 
to be 0.592 ± 0.044 ± 0.018, where the first uncertainty is 
statistical and the second systematic due to uncertainties 
in the correction factors. This measurement is consistent, 
and more precise than 0.79 ± 0.19, the average of previ- 
ous measurements t4| . The measurement for Bj^o ^^ is 
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TABLE VIII. Systematic uncertainties on Fi, \ai\, and <f>i in units of statistical standard deviations (a) for the decay D — J> 
KgK^TT^ . The different contributions are: (I) Blatt-Weisskopf penetration factors; (II) mass and width of resonances; (III) 
quantum correlations; (IV) CLEO III mistag; (V) background models; (VI) sample purities; (VII) acceptance; (VIII) resolution; 
(IX) fitter bias. 









Source {a) 






Total (a) 




I 


II 


III IV V VI VII 


VIII 


IX 




K'-{892)+K- F, 


0.36 


0.33 


0.13 0.08 0.08 0.15 0.16 


0.05 


0.18 


0.59 


K^892fK°s \a,\ 


0.19 


0.24 


0.06 0.09 0.12 0.16 0.14 


0.06 


0.06 


0.42 


K*{892fK''s 0. 


0.87 


0.28 


0.11 0.02 0.10 0.08 0.11 


0.07 


0.05 


0.94 


K''{892fK°s Fi 


0.21 


0.21 


0.05 0.10 0.13 0.19 0.16 


0.05 


0.21 


0.47 


K* {U10)° K^s \a^\ 


4.35 


0.85 


0.23 0.04 0.09 0.07 0.09 


0.01 


0.12 


4.44 


K'{U10fK°s (t>r 


1.13 


0.39 


0.05 0.07 0.17 0.13 0.12 


0.01 


0.13 


1.23 


A'*(1410)°A's Fi 


0.96 


0.35 


0.32 0.04 0.07 0.07 0.06 


0.02 


0.05 


1.08 


K'{1&8Q)+K- |a. 


1.53 


0.41 


0.07 0.08 0.10 0.15 0.14 


0.03 


0.21 


1.62 


K*(IG80)+K' 0, 


1.15 


0.41 


0.03 0.09 0.11 0.17 0.16 


0.02 


0.26 


1.28 


K*{1&8Q)+K- F, 


1.20 


0.32 


0.17 0.00 0.03 0.02 0.12 


0.02 


0.09 


1.26 


K*o{lAm) + K- |a,| 


1.44 


0.32 


0.20 0.08 0.14 0.15 0.21 


0.02 


0.39 


1.57 


i5'o*(1430) + A- 4>^ 


0.49 


0.19 


0.08 0.02 0.08 0.04 0.13 


0.00 


0.05 


0.56 


iCo*(1430)+if- F^ 


1.51 


0.35 


0.15 0.07 0.12 0.13 0.18 


0.03 


0.16 


1.59 


ao(1450)"7r+ \ai\ 


0.41 


0.35 


0.20 0.04 0.08 0.08 0.15 


0.02 


0.08 


0.61 


ao(1450)"7r+ 4>i 


0.42 


0.44 


0.18 0.03 0.04 0.06 0.08 


0.04 


0.06 


0.64 


ao(1450)"7r+ F, 


0.54 


0.24 


0.29 0.02 0.10 0.04 0.22 


0.01 


0.73 


1.01 



TABLE IX. Systematic uncertainties on Fi 
K%K^-K^ . The different contributions are: 
quantum correlations; (IV) CLEO III mistag 
(IX) fitter bias. 



\ai\, and (fn in units of statistical standard deviations (a) for the decay _D" ^• 
(I) Blatt-Weisskopf penetration factors; (II) mass and width of resonances; (III) 
(V) background models; (VI) sample purities; (VII) acceptance; (VIII) resolution; 







Source 


(a) 








Total (a) 




I 


II III IV 


V 


VI 


VII 


VIII 


IX 




A-* (892)- 7^+ F, 


0.23 


0.21 0.20 0.05 


0.03 


0.04 


0.11 


0.00 


0.02 


0.39 


K'{892fK% \a,\ 


0.13 


0.37 0.04 0.03 


0.03 


0.02 


0.05 


0.00 


0.08 


0.41 


K*{892fK% 0, 


0.20 


0.24 0.03 0.01 


0.04 


0.02 


0.12 


0.00 


0.04 


0.34 


K*{892fK% F, 


0.10 


0.23 0.01 0.02 


0.02 


0.04 


0.07 


0.01 


0.07 


0.27 


K5{U30)-K+ \a,\ 


0.33 


0.33 0.10 0.03 


0.02 


0.03 


0.08 


0.00 


0.19 


0.52 


K5{U30)-K+ 4>, 


0.23 


0.43 0.07 0.02 


0.01 


0.02 


0.03 


0.00 


0.13 


0.51 


K^il'i'iG)- K+ F, 


0.31 


0.37 0.09 0.03 


0.02 


0.01 


0.08 


0.00 


0.15 


0.52 


oo(1450) + 7r" \ai\ 


0.58 


0.56 0.08 0.02 


0.01 


0.02 


0.04 


0.01 


0.06 


0.81 


ao(1450) + 7r- 0i 


0.46 


0.58 0.06 0.01 


0.02 


0.01 


0.04 


0.00 


0.04 


0.75 


ao(1450) + 7r- F, 


0.43 


0.24 0.03 0.00 


0.01 


0.01 


0.03 


0.01 


0.17 


0.52 


p(1700)+7r- \ai\ 


1.25 


1.22 0.05 0.00 


0.02 


0.01 


0.07 


0.01 


0.09 


1.75 


p(1700) + 7r- (Pi 


0.36 


0.71 0.09 0.00 


0.01 


0.01 


0.06 


0.00 


0.02 


0.80 


p(1700) + 7r- F 


1.01 


0.68 0.01 0.01 


0.01 


0.01 


0.09 


0.01 


0.02 


1.22 



TABLE X. The results of the Dalitz plot fit for the decay D^ — )• KgK^iv'^ . This is the best model including the neutral 
K''{892) resonance, favored model 1. The amplitudes, phases, and fit fractions are given. The first uncertainty is statistical, 
and the second systematic. 



Amplitude 



Phase 



Fraction (%) 



Component 

D" -^ K*{892)+K- 1.0 (fixed) 0.0 (fixed) 67.6±6.4±3.8 

D" -> A'*(892/'A'g 0.2±0.1±0.0 96.1±13.2±12.4 1.8±1.7±0.8 

D'^ ^ K'luiOf K'^ 8.4±0.7±3.0 -16.7± 7.5 ± 9.2 23.2±3.9±4.2 

D" -^ K'lieSOy K- 15.3±2.3±3.8 148.3±11.1±14.2 22.7±4.4±5.6 

D'^ ^ KSiU30)+K- 1.3±0.3±0.5 -101.0±22.0±12.3 2.5±1.2±2.0 

D" ^a()(1450)~7r+ 0.4±0.2±0.1 ~26.8±42.3±27.2 1.0±0.9±0.9 
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FIG. 5. Fit projections for the D" — >■ i^'^if"?!"'" decay mode using favored model 1, the best model containing the 7^* (892)" 
resonance. The data are shown as points and the total fit projection is overlaid with the solid line. The projection with the 
dotted line is the contribution of background. Plots (a)-(c)((d)-(f)) show the CLEO III (CLEO-c) data for the three pairs of 
invariant mass distributions. 



TABLE XL The results of the Dalitz plot fit for the decay D" — > KgK tt^. This is the model with the lowest x^ per degree 
of freedom, favored model 2. The amplitudes, phases, and fit fractions are given. The first uncertainty is statistical, and the 



second systematic. 



Component 


Amplitude 


Phase (°) 


Fraction (%) 


D" -> K*{S92) + K- 


1.0 (fixed) 


0.0 (fixed) 


62.2±4.8± 3.9 


D° ^ K*{UWfK°s 


10.3±1.0±2.7 


-18.5±5.6± 7.0 


32.0±6.0±19.6 


D° -^ K* (1680)+ K- 


17.1±2.3±2.4 


151.6± 8.3 ± 7.3 


26.1±6.0±24.0 


D° -^ p(1700)-7r+ 


5.5±1.8±1.5 


23.5±12.5±21.5 


4.4±3.0± 3.7 


D" ^KS{U30)+K- 


1.5±0.3±0.7 


-101.5±18.9± 8.7 


3.1±1.2± 3.1 


D° -^ ao(1450)~7r+ 


0.7±0.2±0.2 


-31.7±23.3±16.1 


2.3±1.2± 1.5 



TABLE XII. The results of the Dalitz plot fit for the decay D" — )► KgK^n' . This is the best model including the neutral 
is'* (892)° resonance, suppressed model 1. The amplitudes, phases, and fit fractions are given. The first uncertainty is statistical, 
and the second systematic. 



Component 


Amplitude 


Phase (°) 


Fraction (%) 


D" ^ K*{892)-K+ 


1.0 (fixed) 


0.0 (fixed) 


20.4±2.1±0.8 


D° ^ K* {892 fK°s 


0.4±0.1±0.0 


109.9±10.8±3.7 


3.9±1.5±0.4 


D° ^ K-^{1A?,G)- K+ 


6.2±0.9±0.5 


-114.4± 9.8 ±5.0 


18.4±2.8±1.5 


D° -^ ao(1450) + 7r~ 


3.1±0.5±0.4 


-54.0± 8.6 ±6.4 


15.8±4.1±2.1 


D° -> p(1700) + 7r- 


11.4±2.2±3.9 


90.3± 9.6 ±7.7 


6.2±2.1±2.5 
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FIG. 6. Fit projections for the D" — >■ KgK~Tv'^ decay mode using favored model 2, the best model. The data are shown as 
points and the total fit projection is overlaid with the solid line. The projection with the dotted line is the contribution of 
background. Plots (a)-(c)((d)-(f)) show the CLEO III (CLEO-c) data for the three pairs of invariant mass distributions. 

TABLE XIII. The results of the Dalitz plot fit for the decay D° -^ K%K+-k- . This is the model with the lowest x^ Per degree 
of freedom, suppressed model 2. The amplitudes, phases, and fit fractions are given. The first uncertainty is statistical, and 
the second systematic. 



Component 



Amplitude 



Phase (°) 



Fraction (%) 



DO 



K*{892)-K+ 1.0 (fixed) 0.0 (fixed) 21.1±2.2±3.3 

Kalu30)-K+ 8.0±0.9±0.7 -71.8±13.4± 9.3 31.5±6.2±8.0 

p(1700) + 7r~ 12.7±2.2±5.4 103.7±16.1±11.2 8.0±2.2±3.8 

ao(1450) + 7r- 1.8±0.4±0.3 -30.1±14.7±10.9 5.8±2.3±1.6 

KijUSOfK'^ 4.7±1.5±2.0 --104.1±14.9± 7.5 4.3±2.8±3.4 



also required in the restricted kinematic region of 100 
MeV/c^around the K*{892)^ resonance for the coher- 
ence factor measurement and is here determined to be 
Bk'K = 0.356 ± 0.034 ± 0.007. 



V. COHERENCE FACTOR AND AVERAGE 

STRONG-PHASE DIFFERENCE 

MEASUREMENT 

A. Quantum— Correlated D Decays 

The coherence factor Rj^oj^^ and the average strong- 
phase difference 5j^^ can be determined using double- 
tagged D^D^ pairs from the ^(3770) resonance produced 



in e+e~ collisions at ^/s = 3770 MeV. The two mesons 
are produced in a C-odd eigenstate and therefore the de- 
cays of the two D mesons are quantum-correlated. There 
are three categories of double tags used. The first of these 
are the flavor tags, where the samples of double tags are 
split into two classes: those where the tag kaon and the 
kaon from the signal side have the same sign of charge 
(SS) and those where the two kaons have opposite sign 
of charge (OS). The second category of double tags are 
those where the tag D meson decays to a CP eigenstate, 
which can either be even or odd. The final category of 
double tags used are those where the tag D meson de- 
cays to iir^^7r+7r~. Each category of double tags has a 

different sensitivity to Rj^o ^^ and 5^^ 

The rate for the two neutral D mesons produced 
from -0(3770) where one decays to the K^K^tt^ fi- 
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FIG. 7. Fit projections for the D" -^ KgK^n~ decay mode using suppressed model 1, the best model containing the K*{892)'' 
resonance. The data are shown as points and the total fit projection is overlaid with the solid line. The projection with the 
dotted line is the contribution of background. Plots (a)-(c)((d)-(f)) show the CLEO III (CLEO-c) data for the three pairs of 
invariant mass distributions. 

nal state and the other decays to any final state G, r{KgK^7r^ \G), is given by Q 



T{K'lK^7r^\G)=To 



\-^K°K7rimK'>K^^K7r)-^Gi^) - -^-Eoj^iTnKOK,TnK^)AG{z)\ drrij^oj^drnj^^dz 



fK'iKTT^ 



= ^o[Aj,OK.A^ + Aj^p^A'a - 2Ri,0K.RGAK0K.Aj^^^AGA^cosiS'B - 5^" )]- (16) 

I 

where z describes the multi-body phase space for final is the normalization factor. Specifying that the tag D 

state G, Rg and Sg are the coherence factor and average meson decays to either a flavor tag or a GP tag, Eq. ([T| 

strong-phase difference for the final state G as defined simplifies to the following relations: 
in an analogous way to that given in Eq. (|31), and Fq 



^0.._,^o±^_._.2 ,2 M , ..^g^-^2 o...^^sK^Ji .^^^KOK. 



nK^sKTTlCP^) = ToAUj,^A'cp[l + (r^^ 



- 2\±rr 






m^sKnlF'') - roAW^A%[{4^^y + [r^r - 2r 






^K^'K.^ 



r{K'^KTT\F 



OS\ 



K'-^.K^ F^,2 



^^AUkAf[^ + irn'^'r'^y - 2r^r---"i?^o^,i?^ cos(<5£, - S^^'^")], 



(17) 
(18) 
(19) 



where F, the fiavor tag is either K^tt^, K^n^n^TT , or K^tt^tt'^, the superscript SS or OS denotes the charge 
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FIG. 8. Fit projections for the D — ^ K^K -k^ decay mode using suppressed model 2, the best model. The top plots show the 
projection of the fit on the CLEO III data while the lower plots show the projection onto CLEO-c data. The data are shown 
as points and the total fit projection is overlaid with the solid line. The projection with the dotted line is the contribution of 
background. Plots (a)-(c)((d)-(f)) show the CLEO III (CLEO-c) data for the three pairs of invariant mass distributions. 



relation between the signal and tag kaon, and CP denotes 
a CP eigenstate with eigenvalue A± = ±1. The quantity 
r^ is analogous to Eq. (JU using the DCS and CF am- 
plitudes for the given flavor tag. These rate expressions 
ignore CP violation in D decays, which is well moti- 
vated theoretically and experimentally known to be very 
small [23 . |23| . These rate equations allow for the number 
of expected events in a given data sample to be calculated 
for each category of double tag. The left hand side can 
be replaced with the efficiency-corrected background- 



J 



subtracted signal double-tag yield, S{K^Kti\C), where 

S{KIKtt\G) = N{K%Kn\G)/e{K%KTi\G), (20) 

where N{K^Kn\G) is the background-subtracted yield 
of double tags where one D meson decays to K^Kn and 
the other to final state G, and e{KgKTr\G) is the effi- 
ciency to select that particular double tag. Ratios of the 
rate equations are taken to remove the dependence on 
the normalization factor and Aj^oj^^. This leads to the 

useful observables v^ and k^ defined as 



S{KIKt:\FOS) 



K? + (r^'^'y 



2r 



D' D 



Rf [p cos((5^) + q sin((5'^)] 



l + (rg)2(r^-«-)2_2, 



K%K-!v 



D' D 



RF[pcos{SP) + qsm{SP)] 



(21) 



S{K°Ktt\CP^)B{D ^ Kit) 
^ S{K°KTr\KTrSS)B{D -^ CP) 



1 + (^i;' )' - 2A±r^-^> 



(1 - A±j/)(r^'^)2 + (r^^^"")^ - 2r^-^r'^"^'"[pcos{S'<^) + qsm{5K^)]' 



(22) 



where 



9 = ^A'OK7rSmd^^ , 



(23) 
(24) 



I 

and y is the charm mixing parameter given by y = (r2 
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ri)/2r where r2 and Fi, are the lifetimes of the two D 
meson mass eigenstates and T is the mean of r2 and Fi. 
The branching fractions are included using the following 
relations: B[D^ -^ CP^) = Alp{l T y) and B{D° -^ 
K'TT^) — A^_^+. In the case of the n observable the 
yields and branching fractions are combined so that the 
true value of the k"*" (k~) observable is the same for all 
CP even (odd) tags. 

The other tag group which also provides sensitivity 

to Rxokti and i5^^ ^ is K'^^tt^tt^. The Dalitz plot 
of Kg ^7r+7r~ can be described by two coordinates rn^ 
and m^, where m^{m^) is the invariant mass of the 
Kg j^TT^ {Kg j^TT') pair. The Dalitz plot can be divided 
into bins where the i*'' bin has rri^ > m^ , and the 



symmetric bin in the region where m^ < 



is la- 



beled the 



„-th 



A6d, is defined as A5d = <^_D(m^,m^) — SD{m'^,m'^), 
where 60(171^, m?) is the phase of the D° decay at that 
point. In each bin the amplitude-weighted averages of 
cos(A((5d)) and sin(A((5£))) over each Dalitz plot bin, de- 
noted Ci and Si respectively can be determined. These Ci 
and Si measurements have been performed [14| , together 
with the equivalent parameters for K'^t:'^t:~ decays, cj 
and sj. In addition the flavor tag yields in each bin, Li 
and L^i {L[ and L'_j) for Kg7T~^7r~ (i^^7r+7r~) have also 
been determined. The binning choice used in the current 
analysis is that of the equal regions in A{6ij) as defined 
in [I4i]. 



In quantum-correlated decays of D mesons at ^/s = 
3770 MeV, the yield of the double tag where the tag D 

th 



bin. For each point the phase difference, meson decays to Ks^lt^ it in the i bin is given by 



S{K°gKTT\K''gTT+n-y, = Hs{L, + (rf^'fL^, - 2rf^^{I~L~)[c,p + s,q]), 
S{KlKn\Kln+n'), = Hl^L^ + (r^"^^)^^'^, + 2rf''^^{L[L'_MP + «',9]), 



(25) 
(26) 



r 



where Hs and H^ are normalization factors, and S is 
the efficiency-corrected signal yield in the i**^ bin of the 
Dalitz plot. The inclusion of the K'g ^tttt tags is necessary 
in the measurement as the flavor and CP tags on their 
own provide limited sensitivity to q. 



B. Measurement of observables 

The background-subtracted signal yields of each dou- 
ble tag are given in Sec. Ill Bl To determine the efficiency- 
corrected yields it is necessary to know the efficiency for 
each double tag. The double-tag reconstruction efficien- 
cies are determined from separate simulated signal sam- 
ples for each final state. The resonant decay structure is 
not simulated for any of the D" meson decays because no 
significant differences are observed between the overall ef- 
ficiency determined with flat simulation and the resonant 
decay simulation. In each final state 20,000 events are 
generated, except for the K^ ^tt+tt" tags where 200,000 
events are generated so that the statistical uncertainty 
on the determined efficiency remains small in each bin of 
the Dalitz plot. The determined efficiencies are corrected 
for known small differences between data and simulation 
due to 7r° reconstruction and particle identification effi- 
ciencies. The double tag reconstruction efficiencies range 
from 20% for tt+tt" to 2% for the Kl'q' tag. The vari- 
ation in efficiency from bin to bin in either K'giT'^'n~ or 
K^TT^TT' tags is less than 2%. The effect of the reso- 
nant structure on the determined efficiency is taken into 
account as a systematic uncertainty. 

As the resonant decays of either D meson are not sim- 



ulated the measured efficiencies for the SS flavor tags 
and OS flavor double tags are identical. The system- 
atic uncertainty arising from any difference that does ex- 
ist is considered separately. Hence the v values defined 
in Eq. (|211) are simply the ratios of the background- 
subtracted signal yields given in Table HIl and are 0.66 ± 
0.10, 0.63 ± 0.08, and 0.54 ± 0.06 for u^" , v^^^^ ^ and 
i/^'^'"", respectively. The values of the S{KII£k\CP) and 
S{K%Kt:\Kt:'^'^) yields are fisted in Tables |XlV] and |XVl 
respectively. The values of the S{K^gKTi\K^TnT) yields 
are listed in Tables |XVIl and |XVlIl 

To determine k for each CP tag the branching frac- 
tions are taken from [J]. All subsidiary branching frac- 
tions such as K^ — )■ 7r~^TT~ are included in the calcula- 
tion. The uncertainties on the k value include a statisti- 
cal component from the double-tag yield, and systematic 
uncertainties arising from the uncertainties in the differ- 
ence between Monte Carlo and data, the scale factors 
used to determine background yields, the uncertainties 
of branching fractions used, and the yield of the Kw^^ 
double tag. A number of these uncertainties are corre- 
lated among some or all of the observables. 

The systematic uncertainties due to branch- 
ing fractions, the background scale factors, and 
N{KgKTr\KTr^^) are determined by shifting the de- 
fault value of the branching fraction, scale factor or 
N{KgKn\KTr^^) by one standard deviation and taking 
the shift on the observable as the systematic uncertainty. 
There are systematic uncertainties due to uncorrected 
differences between simulation and data tracking and 
reconstruction efficiencies. Relative uncertainties of 
2% for each tag containing a tt", 0.6% per kaon track. 
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TABLE XIV. Efficiency-corrected yields for double tags with no K^l in final state. The uncertainties take into account the 
statistical uncertainties on the uncorrected yields only. 

l^ag TTK '^ A'^tt^tt" Kfi?^ K%^ A'g77(77) A'"g77(7r7r7r") Kfq' lol^ 

S'(A'5iA'7r|G) 10.7±16.0 16.1±15.1 274.5±89.8 305.0±65.7 293.7±89.9 93.8±36.1 7.5±18.6 52.1±62.6 989.5±110.5 



TABLE XV. Efficiency-corrected signal yields for the CP tags with a K^. The uncertainties take into account the statistical 
uncertainties on the uncorrected yields only. 

Tag i^^» Klu. Kin KW KW^^ 

S{K%Kti\G) 171.6±45.0 58.0±53.0 -1.0±13.1 9.1±34.0 363.2±262.3 



0.3% per charged pion track, and 0.8% for each CP tag 
containing a K^ are applied, determined from studies 
of simulated and real data. The shift in k, due to the 
altered efficiencies is taken as the systematic uncertainty. 
The systematic effect due to non-uniform efficiency over 
the Dalitz plot is treated separately and discussed in 
Sec.lVC] 

To combine the systematic uncertainties an error ma- 
trix is constructed where the diagonal elements are sep- 
arate uncertainties combined in quadrature, and the off- 
diagonal elements are the sum of all ajG}. where the ef- 
fect of the uncertainty is correlated between observables 
j and k. The systematic uncertainties on the k observ- 
ables are given in Tables IXVIIII and IXX| and the cen- 
tral values, statistical, and total systematic uncertainties 
summarized in Tables IXIXI and IXXII The systematic 
uncertainties are smaller than the statistical uncertain- 
ties. The leading systematic uncertainties arise due to 
uncertainties on the branching ratios and the estimated 
background yields. The average values for k+ and k~, 
taking into account systematic uncertainties and corre- 
lations, are 0.56 ± 0.22 and 2.64 ± 0.67, respectively. 
These averages are determined by minimizing the x^ be- 
tween the measured and average k values. The x^/dof is 
12.0/11 which indicates a good overall compatibility of 
the individual measurements. The significant difference 
between k"*" and k^ values indicates the presence of quan- 
tum correlations and a non-zero value for the quantity 
p. In comparison with the incoherent expectation there 
is an excess of CP-odd tags and a reduction in the yield 
of CP-even tags, and hence even the low yields of the 
CP-even tags carry significant information. The error 
matrix for the k observables is given in Appendix B. 

The systematic uncertainty on the v observables is 
very small as there are no contributions from efficiencies 
or branching fractions. The only systematic uncertainty 
considered is that on background yields which partly can- 
cels in the ratio. In the case of the K^ ^tt+tt^ tags the ob- 
servable is simply the efficiency-corrected background- 
subtracted yield. The background yield uncertainty is 
calculated in the same way as for the CP tags and the 
correlations between the uncertainties for the X^tt+tt^ 
tags and the CP tags are determined, and included in 



the error matrix. Global changes to the efficiencies, those 
which are the same for all bins, are not relevant as any 
effect on the efficiency-corrected yield will simply be ab- 
sorbed into the normalization factors in the fit. The 
yields are determined per bin in the Dalitz plot, which 
gives rise to the possibility that detector resolution will 
lead to bin migration. This is not considered here as bin 
migration effects are already taken into account in the 
uncertainties on the measured values of Ci and Si . 



Combination and further systematic 
uncertainties 



To determine the coherence factor and average strong- 
phase difference the information on the observables is 
combined in a x^ fit to determine the best-fit values for 
p and q. To take into account all the uncertainties and 
their correlations the x^ is given by 



X 



ApE-^Ap^, 



(27) 



where Ap is a vector of differences between the expected 
and observed values of the observables and E is the er- 
ror matrix derived from the total uncertainties and their 
correlations. The expected values are computed using 
either Eq. dH]), Eq. (HU, Eq. ^, or Eq. ^ as 
appropriate for a given value of p and q. To calculate 
the expected values a number of external inputs are re- 
quired. These external inputs are free parameters in the 
fit but are Gaussian-constrained to their measured val- 
ues, which takes into account any systematic uncertain- 
ties arising from the measured central values of the exter- 
nal inputs. Correlations between the uncertainties of the 
external parameters are taken into account where known. 
The values of S^"^ , r^^ , x, and y and their uncertainties 
and correlations are taken from |24| . The information 

on Cj- , si is from [lj| , and the coherence information 

on other final states from [7]. The value of r^j^ used 
is determined using the measured value of B^o ^^ and 
the following relation between Bf^o ^^ and the D mixing 



20 



TABLE XVI. Efficiency-corrected signal yields in the K'^-k^-k bins. The uncertainties take into account the statistical 
uncertainties on the uncorrected yields only. 

^iiu i 2 3 4 5 6 7 8 

S{KlK-K\K%-K+-K-)i 55.9±26.0 17.3±15.0 18.0±14.4 30.8±18.1 70.8±28.9 8.8±10.0 20.1±16.8 37.2±20.2 

Bin; -1 -2 -3 -4 -5 -6 -7 -8 

S(K%K-K\Kl-K+-K-)i 42.0±24.1 40.5±22.2 53.6±24.9 35.3±19.4 81.4±29.6 43.2±21.9 5.3±11.1 76.9±29.6 



TABLE XVII. Efficiency-corrected signal yields in the KI-k^-k bins. The uncertainties take into account the statistical 
uncertainties on the uncorrected yields only. 

^liu 1 2 3 4 5 6 7 8 

S{K%K-K\Kl-K+-K-)i 123.6±28.4 28.7±14.5 28.5±13.8 4.0±6.8 20.7±13.2 18.2±12.2 107.0±25.2 106.5±25.5 

Biui -1 -2 -3 -4 -5 -6 -7 -8 

S{K%KTT\KlTT+Ti-)i 128.5±29.4 48.5±17.5 20.1±n.5 9.5±9.2 14.3±12.3 2L3±13.0 53.5±18.7 58.0±19.1 



TABLE XVIII. The systematic uncertainties on the k+ values. 



KK ^^ KItv'' K'loj Kin KU' K'>n\' 



trrtr 



CP tag branching fraction 0.004 0.017 0.095 0.077 0.005 0.060 0.553 

D° -> Ktt branching fraction 0.003 0.012 0.019 0.007 0.001 0.006 0.042 

Background scaling ratio 0.057 0.196 0.062 0.028 0.128 0.090 0.015 

JCtt tag yield 0.028 0.114 0.183 0.064 0.006 0.058 0.408 

CP tag tracking and reconstruction 0.002 0.004 0.029 0.010 0.000 0.002 0.070 

Ktv tracking and reconstruction 0.001 0.006 0.010 0.003 0.000 0.003 0.022 

Total syst. 0.064 0.228 0.218 0.105 0.129 0.122 0.693 



TABLE XIX. The central values k+ and the statistical and systematic uncertainties. 



KK Tvn KW Kluj Klrj Kir)' K'^tt"ti 



Central k value 0.219 0.909 1.456 0.512 -0.049 0.459 3.247 

Statistical uncertainty 0.328 0.855 0.382 0.467 0.665 1.702 1.063 

Total systematic uncertainty 0.064 0.228 0.218 0.105 0.129 0.122 0.693 

Total Stat. + syst. 0.334 0.885 0.440 0.479 0.677 1.706 1.269 



TABLE XX. The systematic uncertainties on the k values. 





AgTT KsU] 


KMll) 


7^g»7(^^^" 


) K'^n' 




Tag branching fraction 


0.145 0.571 


0.773 


0.117 


0.498 


0.499 


D" — >• K-K branching fraction 


0.040 0.049 


0.089 


0.013 


0.049 


0.038 


Background scaling ratio 


0.004 0.015 


0.010 


0.030 


0.008 


0.075 


Kit tag yield 


0.390 0.476 


0.873 


0.131 


0.483 


0.368 


CP tag tracking and reconstruction 


0.067 0.083 


0.063 


0.023 


0.038 


0.057 


K-K tracking and reconstruction 


0.021 0.025 


0.047 


0.007 


0.026 


0.020 


Total syst. 


0.424 0.750 


1.172 


0.180 


0.697 


0.628 
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TABLE XXI. The central values k and the statistical and systematic uncertainties. 





Kgn 


Ks^ 


Ksvi-y-y) 


K'^vi^Tvi:") 


kW 


KU'n^' 


Central k value 


3.102 


3.791 


6.949 


1.043 


3.841 


2.928 


Statistical uncertainty 


0.669 


1.160 


2.676 


2.595 


4.614 


2.114 


Total systematic uncertainty 


0.424 


0.750 


1.172 


0.180 


0.697 


0.628 


Total Stat. + syst. 


0.792 


1.382 


2.921 


2.602 


4.667 


2.205 



parameters 0: 



KSKtt 



^K%K-K 



(rjy'' )2 + {yp 



s K%K-n 



1 + {yp — xq)r 



K°KTr 

D 



(28) 



where the quadratic and higher order terms in y and x 
are neghgible. 

The x^ is minimized to determine the best-fit values 
of p and q. The normalization factors Hs and H^ of 
Eqs. ((25|) and p6|) are also free parameters in the fit. The 
best-fit values are p = 0.719 ± 0.077 and q = 0.105 ± 
0.194. The xVdof of the fit is 53.8/41. The performance 
of the fit is tested using simple simulation. The simple 
simulation generates observable values based on input 

values of Rj^oj^^ and (5^^ , where the observable value 
for each tag is smeared by the experimentally observed 
uncertainties. The studies using the simple simulation 
confirm that the use of tag yields that are consistent with 
zero do not cause bias and that tags with low yields do 
add sensitivity to the overall measurement. 

The efficiencies are determined from Monte Carlo 
where the decay products are generated uniformly over 
phase space. However the efficiency is not uniform over 
the Dalitz plot, and furthermore as the distribution of 
events over the Dalitz plot is different depending on the 
tag due to quantum correlations a further systematic un- 
certainty must be considered. This uncertainty is deter- 
mined by using the isobar models described in Sec. IIIII 
The shifts in the values of p and q are computed using 
Eq. ([3|) and inserting either a fiat efficiency or the CLEO- 
c efficiency model for a variety of isobar models including 
the reported ones. The largest shifts observed between 
the flat and CLEO-c efficiency model are 0.016 for p and 
0.002 for q which are small in comparison to the mea- 
sured uncertainty. Hence these shifts are assigned as the 
systematic uncertainty from non-uniform efficiency. 

As p is positive and Rj^o p^^ is defined to be between 

zero and one, the solution for 5^^ lies in the range 
from -90° to 90' 



and 5 



^s'^" _ 



D 



. This leads to Rk^k^ = 0.73 ± 0.08 

3±15.2)°. Figure IH] shows regions of 

{RK°K-m ^D^ ) parameter space from the fit consis- 
tent with one, two, and three standard deviations from 
the best-fit point. The likelihood is computed as £ = 
g-(x -Xmiii)/2 at a point in parameter space; the fit is re- 
peated at each point with the values of p and q fixed. The 
isobar models also give predictions of the coherence fac- 
tor via Eq. ([3]). The isobar models predict a coherence 
factor of 0.72 (0.79) when the central values of favored 
and suppressed models 1 (2) are used. The predictions 
are consistent with the measured value. 



D. Coherence factor in a restricted region 

Around dominant resonances occurring in both Dalitz 
plots, the coherence factor is expected to be close to one. 
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FIG. 9. The best-fit point for Rj^o j^^ and 5j^^ measured 
over the whole Dalitz plot, and the regions enclosing 1, 2, and 
3 standard deviations from that point. 



as this situation approximates the two-body decay which 
is fully coherent. Higher coherence in principle provides 
increased sensitivity to 7. Therefore the analysis is re- 
peated using only the data where the invariant mass 
of the K''gTi combination lies within 100 MeV/c^ of the 
A'* (892)* mass. The analysis could, in principle, also 
be performed in the region outside this invariant mass 
window; however the data are insufficient for this com- 
plementary study. The efficiencies for each double tag 
are recalculated for the restricted region. The efficiency- 
corrected yields for the CP tags and Ktt^^ tag in this 
restricted region are given in Tables IXXIII and IXXIIII 
The efficiency-corrected yields for Kg ^tt+tt" tags are 
given in Table IXXIVI and Table IXXVl The expected 
yield of the Klrj' tag is very small and is not used. The 
error matrix relating the uncertainties for the k observ- 
ables is given in Appendix B. In the restricted region, 
the values of the ly variables are 0.40±0.08, 0.35±0.06 
and 0.35±0.05 for i^^'^, jy^'^'^^, and ly^"''" , respectively. 

The value of Bk'K that is used in the analysis is the 
one calculated for the restricted region in Sec. IIVI In 
the restricted kinematic region case, a further system- 
atic uncertainty arises due to detector resolution as the 
candidate events may migrate across the defined bound- 
ary in the Dalitz plot. To account for this migration 
simulated signal data samples are generated using the 
resonant models determined in the amplitude analysis. 
For the CP tags and the K'g j^7r"'"7r~ tags the quantum 
correlation between the tag and signal side is emulated. 
Candidate events can migrate into the restricted region, 
or they can migrate out. There is no systematic uncer- 



23 



TABLE XXII. Efficiency-corrected yields in the restricted region for double tags with no /•i'£ in final state. 



KK 



Tag 

S{K"sK-k\G) 6.0±f4.2 2.0±II.I 43.7±37.8 233.1±57.9 f99.9±68.7 47.8±28.2 



K'-^nill) KM^nn^') 



'kW 



Kit SS 



fO.I±I8.I 83.0±87.5 288.9±48.7 



TABLE XXIII. Efficiency-corrected signal yields in the re- 
stricted region for the CP tags with a K^. 



Tag 



liTn Kl^^ 



S{K'-i;KTv\G) 116.2±37.1 25.9±38.9 -4.7±8.5 186.4±194.0" 



tainty if the migration in both directions is equal, but 
the presence of the resonance means that a larger frac- 
tion of events will migrate out of the restricted region 
than migrate in. The net migration factor is determined 
from the simulated data. It is different for different types 
of tags due to the quantum correlations which alter the 
distribution of the K^Kir decay on the Dalitz plot. The 
net migration factor is treated in the same way as a sys- 
tematic uncertainty on the efficiency in the analysis. The 
net migration factor, that is the net loss of events inside 
the region, varies between tag categories and is of order 
1% of the yield within the region. 

Tables IXXVlHXXIXI give the results for the observable 
K and its uncertainties for the K*{892)^ bin. The average 
values are k+ = 0.44 ± 0.30 and k" = 3.18 ± 1.08 with a 
X^/dof of 6.5/10. The uncertainties due to non-uniform 
efficiency are 0.002 and 0.001 for p and q, respectively. 
They are reduced in comparison to the full Dalitz plot 
measurement as the efficiency is more uniform over the 
restricted region. 

The fitted values of p and q are p — 0.893 ± 0.089 and 
q = 0.446 ± 0.304, with xVdof given by 48.7/40. This 
leads to a measurement of -Ri<-. if = 1.00±0.16 and 5|^ ^= 
(26.5 ± 15.8)°. The regions consistent with one, two, 
and three standard deviations from the best-fit point are 
shown in Fig. [TU] The coherence in the region of the 
/'ir*(892)^ is larger than the coherence over the full Dalitz 
plot and is consistent with one. The prediction from 
the isobar models for the restricted region is 0.93 (0.94) 
when favored and suppressed models 1 (2) are used. The 
predictions are consistent with the measured value. 



E. Impact on 7 measurements 

It is instructive to estimate the impact of the knowl- 
edge of the coherence factor and average strong-phase 
difference on a measurement of 7 performed at LHCb or 
a future high luminosity e^e~ experiment. Simple sim- 
ulated data are generated assuming 15000 B~^ — 7> DK^ 
events with D decaying to KgK^n^ , and 7 — 70°, Sb — 
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FIG. 10. The best-fit point for Rk*k and 5^ ^ which are 
measured over the restricted region, and the regions enclosing 
1, 2, and 3 standard deviations from that point. 



ues for 7, Sb, and rg are consistent with current knowl- 
edge [1| . The data samples are fit with only the parame- 
ter 7 free. The width of the best-fit 7 distribution is 10°. 
The restricted region will have approximately 60% of the 
statistics of the full Dalitz plot. A similar study in the 
restricted region assuming 9000 events and the values of 

rK'K 



as measured here for the full Dalitz plot. The chosen val- 



and Rk'K and 5^ ^ as measured in the restricted 
region has a best-fit 7 distribution of 5°. The better 
sensitivity, despite the smaller sample, is due to a com- 
bination of increased coherence, a lower value of r^*^ 

in comparison to r^^ , and a different value of the av- 
erage strong-phase difference. The uncertainties on the 
coherence factor and average strong-phase difference are 
incorporated by making them free parameters in the fit 
and applying Gaussian constraints based on their mea- 
sured uncertainties. In this case the width of the best-fit 
7 distribution increases to 14° and 7° for the full Dalitz 
plot and restricted regions, respectively. The uncertain- 
ties Rxoktt, ^d^ '^ ■• Rk'K, and 5^ ^ are dominated by 
the statistical uncertainties and hence the precision of 
any measurement of 7 could be increased with a more 
accurate measurement of the coherence factors and aver- 
age strong-phase differences. The inclusion of this decay 
channel in a global fit to 7 with other channels will im- 
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TABLE XXIV. Efficiency-corrected signal yields in the restricted region in the K%-k^tx bins. 

^h^ 1 2 3 4 5 6 7 8 

S{KlK-K\Kl-K+-K-)i 40.3±36.8 48.3±35.8 19.5±23.6 49.3±35.8 128.2±60.4 25.1±27.5 54.6±44.3 18.8±23.8 

Bin -1 -2 -3 -4 -5 -6 -7 -8 

S{KlKTT\KlTT+-K-)i 65.7±43.7 45.3±38.4 59.0±40.0 23.2±25.4 108.6±50.3 54.8±39.8 16.3±21.8 115.7±59.2 



TABLE XXV. Efficiency-corrected signal yields in the restricted region in the A'£7r+7r bins. 

^m i 2 3 4 5 6 7 8 

S{KlK-K\Kl-K+TT-)i 20L0±54.9 23.3±22.2 60.7±29.2 -9.4±14.0 25.6±2L3 -1.2±13.0 156.2±45.8 176.2±48.5 

Bin -1 -2 -3 -4 -5 -6 -7 -8 

S{KlK-K\Kl-K+-K-)i 242.5±59.6 79.2±33.0 26.5±19.4 2.4±12.6 29.0±22.5 12.9±19.8 86.4±34.7 65.2±29.6 



TABLE XXVL Systematic uncertainties on the k^ values determined in the restricted region. 



KK ^+^- KW' KIlo KU KW*^^ 



Tag branching fraction 0.004 0.004 0.109 0.058 0.045 0.149 

D° -^ K-K branching fraction 0.003 0.002 0.021 0.005 0.005 0.011 

Background scaling ratio 0.091 0.230 0.062 0.032 0.126 0.017 

K-K tag yield 0.042 0.039 0.338 0.078 0.082 0.177 

CP tag tracking and reconstruction 0.002 0.001 0.033 0.008 0.000 0.019 

K-K tracking and reconstruction 0.001 0.001 0.011 0.003 0.003 0.006 

Bin migration 0.000 0.000 0.003 0.001 0.001 0.001 

Total syst. 0.100 0.233 0.363 0.103 0.157 0.233 



TABLE XXVIL The central values of k^ determined in the restricted region and the statistical and systematic uncertainties. 



KK K+K- kW' Klu Kin J^s^"^" 



Central n value 0.205 0.192 1.665 0.386 -0.402 0.873 

Statistical uncertainty 0.491 1.059 0.532 0.579 0.730 0.755 

Total syst. 0.100 0.233 0.363 0.103 0.157 0.233 

Total Stat. + syst. 0.501 1.085 0.644 0.588 0.746 0.790 



TABLE XXVIII. Systematic uncertainties on the k values determined in the restricted region. 

Tag branching fraction 0.188 0.656 0.665 0.266 1.339 0.432 

D° -> Kk branching fraction 0.051 0.056 0.077 0.031 0.133 0.033 

Background scaling ratio 0.004 0.008 0.016 0.045 0.000 0.054 

Kk tag yield 0.812 0.884 1.212 0.483 2.094 0.515 

CP tag tracking and reconstruction 0.086 0.095 0.054 0.052 0.102 0.050 

Kk tag tracking and reconstruction 0.027 0.029 0.040 0.016 0.069 0.017 

Bin migration 0.006 0.005 0.007 0.003 0.011 0.003 

Total syst. 0.840 1.107 1.386 0.557 2.493 0.677 
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TABLE XXIX. The central values of k determined in the restricted region and the statistical and systematic uncertainties. 







K'iu 


7^^,7(77) K'^i^i^nn") 


Kh' 




Central k value 


4.003 


4.357 


5.975 


2.382 


10.324 


2.537 


Statistical uncertainty 


0.994 


1.498 


3.523 


4.269 


10.881 


2.640 


Total syst. 


0.840 


1.107 


1.386 


0.557 


2.493 


0.677 


Total Stat. + syst. 


1.301 


1.862 


3.786 


4.305 


11.163 


2.726 
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prove the overall sensitivity to this parameter. 
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Appendix A: The interference fit fractions 

This appendix gives the full interference fit fractions 
as defined in Eq. (fT4)) between each resonance for each 
model in Tables IXXXtiXXXml 



Appendix B: The covariance matrices 



Tables IXXXIVI and IXXXVI give the the covariance 
matrices for the uncertainties between the k measure- 
ments from the different tags. Entries below the diagonal 
are not shown. 
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TABLE XXX. 
as a percent. 



The interference fit fractions between every pair of resonances in favored model 1. The fractions are provided 



ii'*(892)" 



K*{14W)" 



K*{1680y 



7^,^1430)" 



ao(1450) 



T" 



if* (892)- 
K*{892;f 
ii"*(1410)° 
K*{168(l)- 

K*o{im)- 



-0.23±0.65±0.57 -12.29±1.23±3.42 13.62±3.32±2.17 3.38±0.88±0.63 
-2.81±1.55±0.79 0.23±0.27±0.75 0.57±0.30±0.24 
-21.18±2.73±8.66 1.15±0.61±1.39 
-6.70±1.79±2.6: 



2.56±2.09±1.71 
0.26±0.30±0.12 
0.78±3.19±1.99 
0.55±4.08±2.00 
1.37±1.17±0.75 



TABLE XXXL The interference fit fractions between every pair of resonances in favored model 2. The fractions are provided 
as a percent. 



A'*(1410) 



TT 



K*(1680)- 



p(1700) 



"F" 



7^5(1430)- 



a^^{U50y 



K*{892)- 
J^*(1410)° 
A-* (1680)- 
p(1700)+ 



^14.04±1.42±2.02 13.05±2.61±4.02 3.60±1.52±3.36 
-27.47±5.29±28.07 -7.33±3.45±8.73 
-4.33±2.76±7.87 



3.58±0.70±0.92 
1.55±0.74±2.30 
-8.06±1.93±7.16 
1.24±0.55±1.57 



4.03±1.60±1.52 

1.82±3.41±1.72 

-0.30±4.16±1.9: 

0.02±0.01±0.01 

2.49±1.34±0.86 



TABLE XXXII. The interference fit fractions between every pair of resonances in suppressed model 1. 
provided as a percent. 



The fractions are 



K*{892y' 



KSjusoy 



ao(1450)^ 



p(1700)^ 



ii'*(892)- -0.74±0.46±0.10 4.91±0.57±0.53 7.01±0.88±0.87 5.74±0.91±1.16 

if* (892)" 1.34±0.54±0.33 -0.50±0.58±0.40 0.72±0.47±0.29 

is:o(1430)- 17.39±1.52±1.69 -0.49±0.75±0.41 

ao(1450)+ -0.01±0.01±0.01 



TABLE XXXIII. The interference fit fractions between every pair of resonances in suppressed model 2. The fractions are 
provided as a percent. 





^-1^(1430)- p(1700)+ 


ao(1450)+ iC2*(1430)" 


A'*(892)- 


5.70±0.55±0.31 6.41±0.89±1.41 


3.62±0.95±0.64 1.75±0.81±3.32 


7^0* (1430)- 


2.07±0.85±0.80 


15.25±2.21±0.91 -5.79±5.60±9.15 


p(1700)+ 




0.00±0.01±0.01 4.00±1.23±1.71 


ao(1450)+ 




-3.71±1.64±2.26 



TABLE XXXIV. The covariance matrix for the k observables determined using data from the full Dalitz plot. 





^^^ 


k"" 


.c^-" 


K'^^'" 


^K^rii-n) 


n''^^' 


,x,.%'' 


.^-" 


^Ks^ 


^KsTi{-n) 


^Ksrii-iT-iT-iT") 


^Ks^' 


.^-"-" 


K^^ 


0.11 


0.01 


0.01 


0.00 


0.01 


0.01 


0.01 


0.01 


0.01 


0.02 


0.01 


0.01 


0.01 


k"" 




0.78 


0.03 


0.01 


0.03 


0.02 


0.05 


0.05 


0.06 


0.10 


0.02 


0.06 


0.06 


«^--" 






0.19 


0.01 


0.01 


0.02 


0.08 


0.07 


0.09 


0.16 


0.03 


0.09 


0.07 


«^^" 








0.23 


0.00 


0.01 


0.03 


0.03 


0.08 


0.06 


0.01 


0.03 


0.03 


I^Kr,vill) 










0.46 


0.01 


0.00 


0.00 


0.00 


0.01 


0.01 


0.00 


0.01 


AC^^"' 












2.91 


0.03 


0.02 


0.03 


0.05 


0.01 


0.06 


0.03 


^KsT^vr" 














1.61 


0.16 


0.20 


0.36 


0.05 


0.20 


0.43 


^Ks.° 
















0.63 


0.19 


0.35 


0.05 


0.19 


0.15 


I^Ks^ 


















1.91 


0.42 


0.06 


0.23 


0.18 


^Ksvhl) 




















8.53 


0.21 


0.43 


0.33 


ifs)7(vr7r7r") 






















6.77 


0.06 


0.05 


^Ksv' 
























21.78 


0.18 


«^-"-" 


























4.86 



to 

00 



TABLE XXXV. The covariance matrix for the k observables determined using data from the restricted region. 





^^^ 


k"" 


K^--" 


K'^^'" 


^Kj^ 17(77) 




.--" 


«xs- 


^Ks,,(77) 


^Ksvi-^TTT,") 


K^S'7' 


^K,.'V.' 


^^^ 


0.25 


0.02 


0.02 


0.01 


0.01 


0.01 


0.03 


0.04 


0.05 


0.02 


0.09 


0.03 


k"" 




1.18 


0.03 


0.01 


0.03 


0.01 


0.03 


0.04 


0.05 


0.03 


0.08 


0.03 


.^-^ 






0.41 


0.03 


0.04 


0.06 


0.28 


0.30 


0.41 


0.17 


0.71 


0.18 


K^^" 








0.35 


0.01 


0.01 


0.06 


0.11 


0.10 


0.04 


0.16 


0.04 


^KLnhi) 










0.56 


0.02 


0.07 


0.07 


0.13 


0.06 


0.17 


0.05 


^Ks.o.o 












0.62 


0.14 


0.16 


0.22 


0.09 


0.37 


0.16 


^Ks.^ 














1.69 


0.72 


0.99 


0.39 


1.71 


0.42 


^Ksu. 
















3.47 


1.08 


0.43 


1.86 


0.46 


rv 


















14.33 


0.77 


2.55 


0.63 


^KsviTT^T") 




















18.53 


1.02 


0.25 


^Ksv' 






















124.62 


1.08 


^K^.o.o 
























7.43 



to 

CO 



